void  verletLinkCellList(caja *bx,double Dt)
{
	double posNew[3], dr;
	int index =0;
	int mc[3], c;
	//Reset the head of the link-cell list
	for (int i=0;i<bx->ncell3;i++)
	{
		bx->head[i] = -1;
	}

	for(int i=0;i<bx->npart;i++)
	{
		for (int j=0;j<3;j++)
		{
			//posNew[j]=(bx->part[i].pos[j])+(bx->part[i].vel[j])*Dt+(bx->part[i].force[j])*Dt*Dt/2.0;
			//dr=posNew[j]-(bx->part[i].posOld[j]);
			//dr-=int(dr/(0.5*(bx->lbox[j])))*(bx->lbox[j]);

			bx->part[i].pos[j]+=0.1*(double(rand())/RAND_MAX-0.5);

			//bx->part[i].vel[j]=(dr)/(2.0*Dt);
			//bx->part[i].posOld[j]=bx->part[i].pos[j];
		}
		bx->boundary(i);

		///Creamos la link-cell list
		for(int j=0;j<3;j++)
		{
			mc[j] = int((bx->part[i].pos[j]+0.5*bx->lbox[0])/bx->rcut);
		}
		c = int(mc[0]*bx->ncell2+mc[1]*bx->ncell+mc[2]);
		///Link to the previous occupant (or EMPTY if you're the 1st)
		bx->link[i] = bx->head[c];
		///The last one goes to the header"""
		bx->head[c] = i;
	}
}
